Supersymmetry theory of microphase separation in homopolymer— oligomer mixtures 
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Mesoscopic structure of the periodically alternating layers of stretched homopolymer chains sur- 
rounded by perpendicularly oriented oligomeric tails is studied for the systems with both strong 
(ionic) and weak (hydrogen) interactions. We focus on the consideration of the distribution of 
oligomers along the homopolymer chains that is described by the effective equation of motion with 
the segment number playing the role of imaginary time. Supersymmetry technique is developed to 
consider associative hydrogen bonding, self-action effects, inhomogeneity and temperature fluctu- 
ations in the oligomer distribution. Making use of the self-consistent approach allows to explain 
experimentally observed temperature dependence of the structure period and the order-disorder 
transition temperature and period as functions of the oligomeric fraction for systems with different 
strength of bonding. A whole set of parameters of the model used is found for strong, intermediate 
and weak coupled systems being P4VP-(DBSA)^, P4VP-(Zn(DBS)2)^ and P4VP-(PDP)^, respec- 
tively. A passage from the formers to the latters shows to cause crucial decrease of the magnitude 
of both parameters of hydrogen bonding and self-action, as well as the order-disorder transition 
temperature. 

PACS numbers: 36.20.-r, 64.60Cn, ll.SOPb 
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I. INTRODUCTION 

Surfactant-induced mesomorphic structures based on 
the association between flexible homopolymers and 
head-functionalized oligomers represent a new class of 
supramolccular materials. They exhibit a rich phase be- 
havior to be an object of investigations that have at- 
tracted, during past decade, considerable attention of 
both experimentalists - Q and theoreticians 0,0- 
Microphase separation is the principal property of such 
systems which results in the formation of ordered meso- 
scopic structures due to the association between the head 
group of the oligomer and corresponding groups of the 
homopolymer, on the one hand, and unfavorable polar- 
nonpolar interactions between the non-polar tail of the 
surfactant molecules and the rest of the system, on the 
other one. 

The homopolymer-oligomer systems involve two main 
classes that are relevant to strong ionic bonds and weak 
hydrogen ones. Unlike to conventional copolymers where 
repulsive blocks are bonded together by covalent bonds, 
there are various temporary physical interactions which 
play a crucial role in the formation of ordered mcsophases 
in such systems. In the ionic bonding systems the degree 
of association is relatively high, so that the polymer chain 
resembles a comb copolymer with regularly alternating 
oligomer side chains. At the same time, for the systems 
with temperature-dependent hydrogen bonds the incom- 
patibility must not be so strong to induce separation on 
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a macroscopic level. Here, the microphase separation re- 
sults in the periodic alternation of the layers of stretched 
homopolymer chains surrounded by perpendicularly ori- 
ented oligomer tails (see Figure QJ. Similarly to the con- 
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FIG. 1: Schematic picture of the homopolymer-oligomer mi- 
crophase separated structure for x = 1/3. 

ventional copolymer systems, a rich variety of morpholo- 
gies (lamellar, cylindrical, spherical etc.) is shown to be 
possible n|. However, for the sake of simplicity we will 
restrict ourselves with considering lamellar morphology 
only. 

An example of the ionic bonding systems is represented 
by the mixture P4VP-(DBSA)2: of the homopolymer be- 
ing atactic poly(4-vinyl pyridine) (P4VP) and the sur- 
factant as dodecyl benzene sulfonic acid (DBS A). Here, 
owing to the very strong interaction, the temperature do- 
main of microphase separation is not bounded from above 
by association effects 0, Q . The peculiarity of the sys- 
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terns of this type, being polyelectrolyte-surfactant com- 
plexes, is that the long space structure period is an in- 
creasing function of the number oligomer /monomer ratio 
X (the number of DBSA-groups per one pyridine ring). 
More complicated behavior is inherent in the hydrogen 
bonded systems which were considered to study the oppo- 
site weak-bonding limit ~ 0- Here, the weak interac- 
tion causes an order-disorder transition to homogeneous 
high-temperature state. An example of these systems 
is given by the mixture P4VP-(PDP)2. of the same ho- 
mopolymer P4VP with 3-pentadecyl phenol (PDP) being 
the oligomer. In this case, unlike to the ionic bond- 
ing systems, the long space period decreases with x- 
increase. An intermediate behavior exhibit the system 
P4VP-(Zn(DBS)2)a; with the oligomer being zinc dodecyl 
benzene sulfonate Zn(DBS)2 that forms transition metal 
coordination complexes with the monomers of P4VP . 
Ionic bond weakening due to the absence of covalently 
bound charges along the homopolymer chain leads here 
to a non-monotonic form of the x-dependence of the long 
space period. 

Principally important for our consideration is decreas- 
ing form of the temperature dependence of the long 
space period for all above systems - However, 
such character of the dependence appears in hydrogen 
bonded systems only within a finite temperature interval 
bounded by the glass transition temperature Tg from be- 
low and order-disorder transition temperature Tc from 
above [^Q. Here, an increase of the oligomer/monomer 
ratio X leads to a non-monotonic behavior of the tem- 
perature Tc with a maximum near the point x = 0.85, 
deviation from which narrows the temperature domain 
Tg-h-Tc- This domain is the region of our interest where a 
purely microphase separated structure is possible. Below 
the glass transition temperature Tg the crystallization of 
the oligomer chains occurs that causes a reduction of the 
overall volume of the system and a sudden decrease of 
the long space period 

Microphase separation phenomenon had been exten- 
sively studied in the past two decades for a variety of 
poly mer systems including random heteropolymers [S| - 
[lO| . Theoretical studies of the homopolymer-oHgomer 
mixtures, being the systems of associating polymers, were 
proposed by Tanaka et al. and Dormidontova et al. 
Q within the random phase approximation introduced 
by Leibler Here, the total free energy 

Ftot^Fho + Fhb (1) 

is written as a sum of two terms, Fho related to the non- 
associated homopolymer-oligomer mixture and Fhb at- 
tributed to the hydrogen bonding. Then, making use of 
minimization principle with respect to the dependence of 
the free energy Ftot on the average fraction of hydrogen 
bonds X present in the system, permits to find the tem- 
perature dependence X{T) and to study possible forms 
of phase diagrams for both macrophase and microphase 
separations. It turned out that this approach gives the 
real dependence of the long space period L of the ordered 



structure on the oligomer/monomer ratio in the system, 
however, as the fraction of hydrogen bonds monotoni- 
cally decreases with temperature, the increasing temper- 
ature dependence of L(T) obtained is in contradiction to 
the experimental data [J. This inconsistency is caused 
obviously by the roughness of the random phase approx- 
imation used for description of the hydrogen bonding. 

To avoid this limitation, our approach is based 
on the above mentioned analogy between associat- 
ing homopolymer-oligomer mixtures and random comb 
copolymers taking into account the varying number of 
oligomers attached to the main chain stochastically. Such 
a system can be analysed in terms of the random walk 
statistics to apply the field theoretical scheme JJJj for the 
development of the microscopic theory. The corner stone 
of our approach lies in the assumption that the alter- 
nation of the homopolymer associative groups with and 
without oligomers attached is like the alternation of the 
segments of different types along the chains of a random 
heteropolymer to be represented as a stochastic variation 
of the Ising spin, for which the role of imaginary time is 
played by the number of chain segment ri 12j - [3 | . 

Along this line, the problem under consideration is di- 
vided into two parts, the first of which is reduced to the 
determination of the relation between the long space pe- 
riod L and the average fraction of hydrogen bonds X, 
whereas the second one is focused on the determina- 
tion of the frequency u — 2'nX in the distribution of 
the oligomer heads along the homopolymer chain. The 
first part of the problem was studied on base of the sim- 
plest model that is reduced to the treatment of the 
dependence Fho{L) given by the first term of the free en- 
ergy (Q. Corresponding consideration developed within 
the framework of the strong segregation limit derives to 
generic relation (|A8I) for the dependence L(ijj) (see Ap- 
pendix A). In this paper, we focus on the second problem 
to be related to the definition of an optimal frequency lo 
that minimizes the second term of the free energy ^ 
within the framework of the weak segregation limit. 

The formal basis of our treatment is the field theoret- 
ical scheme of stochastic systems with using the super- 
symmetry field 11]. Conformably to the polymers, this 
theory was proposed in J l5| | an d developed for the random 
copolymers in Re fs. |l2l| - |14(. Our approach is based on 
the Martin-Siggia-Rose method of the generating func- 
tional 01 . Power and generality of the supersymme- 
try field scheme was demonstrated for the Sherrington- 
Kirkpatrick model for which it is identical to the replica 
approach TT]. The formal base of the supersymmetry 
is a nilpotent quantity which represents a square root 
of 0. In this sense, the superfield is similar to the com- 
plex field, in which the imaginary unit, being square root 
of —1, is used instead of the anticommuting nilpotent 
quantity being the Grassmann variable. By definition, 
the supersymmetry field combines commuting the boson 
and anticommuting fermion components into the unified 
mathematical construction representing a vector in the 
supersymmetry space. Choice of the optimal basis of the 
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supersymmetry correlator yields in optimal way the ad- 
vanced/retarded Green functions and the structure factor 
to obtain microscopic expression for the frequency oj. 

The paper is organized in the following manner. Sec- 
tion II contains initial relations of the field scheme used 
to write the system Lagrangian. It involves the effec- 
tive potential energy whose quadratic term describes hy- 
drogen bonding between the oligomers and the associa- 
tive groups of the homopolymer chains, whereas cubic 
and biquadratic terms relate to the self-action effects. 
The principal peculiarity of our approach lies in account- 
ing of the inhomogeneity in the distribution of oligomers 
along the homopolymer chains. This accounting is caused 
by the introduction of the effective kinetic energy whose 
density is proportional to the square of the derivative of 
oligomers distribution over segment numbers n. Due to 
the temperature dependence of the hydrogen coupling, 
related effective mass is a fluctuating parameter whose 
averaging, along the Hubbard-Stratonovich procedure, 
arrives at the biquadratic term with respect to the time 
derivative. According to the calculations given in Section 
IV, just this term, being considered within the mean- 
field approach, causes decaying character of the tempera- 
ture dependence L{T) of the structure period. Complica- 
tion of the problem arising from the determination of the 
proper frequency lj is caused by an essential non-linearity 
and coupling the advanced/retarded Green functions and 
the structure factor. Hence, it is methodically convenient 
to use the supersymmetry technique that enables to ob- 
tain in the simplest way explicit expressions for above 
functions in the long-range limit (see Section III). Di- 
vergency condition of the Green function permits to find 
the proper frequency oj with accounting self-action effects 
within supersymmetry perturbation theory. A compari- 
son of the dependencies obtained with experimental data 
given in Section V shows that the scheme developed al- 
lows to present in a self-consistent manner main peculiar- 
ities of the microphase separation in the homopolymer- 
oligomer systems with associative coupling. 



II. GENERIC FORMALISM 

The problem under consideration is addressed to the 
definition of the effective law of motion c[n) that deter- 
mines a sequence of oligomer alternation along the ho- 
mopolymer chain by means of specifying the occupation 
number being c(n) = 1 if oligomer is attached to the seg- 
ment n, and c(n) = otherwise. When the index of the 
homopolymer chain — s- oo, the argument n may be 
considered as a continuum one, and we are ventured to 
start with Euler equation 



where dot denotes derivative with respect to the segment 
number n, action S and dissipative functional R take the 
usual forms 

N N 

S{c{n)}= [ L{c{n),c{n))dn, ! (c(n))^dn (3) 



being defined by the Lagrangian L{c.[n),c{nj) and the 
damping coefficient 0, respectively. The total action 
5 = AT — n is determined by a "kinetic" contribution K 
of inhomogeneity in the oligomers distribution and "po- 
tential" component 11 = Vq -|- caused by the interaction 
between homopolymer and oligomers 



N 



(4) 



and self-action contribution 

N 

V = T j v{c{n))dn, 





_ M 3 



4!^ 



(5) 



Here, T is temperature measured in energy units, fac- 
tor r determines the strength of the hydrogen bonding, 
multipliers /i, A are self-action parameters. 

In comparison with the above standard approach, the 
determination of the contribution of inhomogeneity along 
polymeric chain is much more delicate problem. Indeed, 
the bare magnitude can be written in the form of the 
usual kinetic action 



N 



K ^T- 



dc 
dn 



dn 



(6) 



where an effective mass m appears as a temperature 
fluctuating parameter with mean value m and variance 
(m — m)^ = (T^ (bar denotes the average, as usually). 
Then, after averaging exponent exp(— AT/T) over the 
Gaussian distribution of the bare mass m, we obtain the 
effective kinetic action in the following form: 



N 

— ~ — Tfl f 9 

IC^K + K; K = T— {c{n)) dn, (7) 



N N 



K = -T- 



{c(n)f{c(n')fdndn'. (8) 



5S_ d_5S_ _ 6R 
5c dn 5c 5c 



(2) 



As a result, total action takes the final form 
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N ^ AT Af AT 

S^T^^J (c(n))'dn-ry J J (c(n))'(c(n'))'dndn' - J (c(n))'dn-y (9) 





where self-action potential V is given by Eqs. ||SJ). Re- 
spectively, Euler equation (O arrives at the equation of 
effective motion 



(m - Aj 



e 



A ) C + — C + TC ~ —v' 



(10) 



where one notices 



Sc{n) 



dv r cr^ 



N 



(c(n'))^dn'. (11) 



By introducing the effective mass m, characteristic num- 
ber of correlating segments ric and ^-correlated source 
C{n) in accordance with definitions 



m = m-A, nc=f, (12) 
(C(n))=0, (C(n)C(n')) =<5(n-n'), (13) 



one obtains Langevin equation of inertial type 

rhii + ricC — —{tc + v') + C.. (14) 

Making use of the field scheme TTl allows to express 
the noise C in terms of an additional degree of freedom 
p being the momentum conjugated to the effective coor- 
dinate c. Following this line, one has to introduce the 
generating functional 



Z{c{n)} ^ {mc + UcC + tc + v' - (} dot 



(15) 

being the average over the noise C(j^) where 5-function 
accounts for the equation of motion (|14|) , and the deter- 
minant is Jacobian of transition from C('^) to c(n) that 
is equal 0/T = Uc. Then, making use of the fmictional 
Laplace representation of (5-function over a ghost field 
p{n) and averaging Eg. ljlSII over Gaussian distribution 
being related to Eas. H13|) . we derive to the standard form 

M 



Z{c{n)} = J exp [-S{c{n),p{n)}]Bp{n), 
S{c{n),p{n)} = C{c{n),p{n))dn 



where effective Lagrangian is introduced 

£ = {mc + TicC + TC + v'){n;^^p) -~ ^{ric^pf- 



(16) 
(17) 

(18) 



According to Euler equations 

dC_d_dC ^dC_^ 

dx dn dx d-n? dx ' ' 

effective motion in the phase space is determined by the 
system 



mc -I- UcC = — (tc + v') + [n^ ^p), 
mp — ricP = — (t + v") p. 



(20) 
(21) 



A comparison of the first of these equations with Eq. H14() 
shows that the conjugated momentum p appears as the 
most probable value of the renormalized noise amplitude 

n-cC 



t pr 

M 



III. SUPERSYMMETRY REPRESENTATION 
OF CORRELATION IN OLIGOMER 
DISTRIBUTION 

Eas. (|20(l . (|21|l represents a system of nonlinear equa- 
tions whose solution demands a use of the perturbation 
theory with respect to the self-action parameters A, fi 
and of the self-consistency procedure to determine an 
effective mass rh{c{n)}. However, because we are in- 
terested in the knowledge not of laws of motion c(n) 
and p{n), but only of the frequency of oligomer alter- 
nation along the homopolymer chain, it is appropriate 
to restrict ourself to an investigation of the correspond- 
ing correlators. The latters reduce to autocorrelator 
S{n) — (c^ln)), and retarded and advanced Green func- 
tions, G-{n) = (c(n)p(O)) and G'+(n) = (p(n)c(0)) re- 
spectively. As shows the consideration given in 12J - [lj|, 
it is convenient to represent these correlators as the com- 
ponents of the supercorrelator 



C(z) = (02(z)), z = {n,d} 
of pseudovectors of the phase space 

= c + {n'^^p)'d 



(22) 



(23) 



being built by making use of nilpotent variable which 
satisfies the conditions: 

^ 0, M' = J d§^0, j §dd = l. (24) 

Along this line, the supercorrelator H22|l appears as a 
pseudovector 



C = G+A + G_B + 5T 



(25) 
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TABLE I: 



l\r 


T 


A 


B 


T 





T 





A 





A 





B 


T 





B 



spanned on set of the orts 

A(i?, d') = t9, B(i?, -d') = d', T(i9, = 1. (26) 

Introducing the functional product of some vectors X, 
Y, Z in such a space: 



(27) 



it is easy to see that orts H26f) are noncommutative to 
obey to the multiphcation rules given in Table I. As a 
result, making use of the supercorrelator presents a 
big advantage in analytical calculations. 

Under suppression of the inhomogeneity fluctuations 
along the homopolymer chain (cr = 0), the action H17() 



with the Lagrangian H18|l written within the harmonic 
approximation (u(c) = const) takes the canonical form 



Mm} = 2 / <t>iz)Lizmz)dz 



(28) 



with the linear operator 



L{z) ^ T{n) + D{z); T{n) = T + fhdl, (29) 



dd I dn 



This operator defines the bare supercorrelator according 
to the relation 

C^'^\z) = L-\z)5{^,{>'), 5{-dJ')^d + d'. (31) 

Taking into account condition = n^d^ '^^^ ob- 
tains 



^(0) ^ [T{n)-D]d{{},r) 
T^{n)-nldl 



(32) 



Using Fourier transformation over the frequency we 
arrive to the expression 



Ci^)^ ^ + [r{-)--^nM^+[r{^) + -^nM^' r(.)^r-m.^. (33) 

Then, taking into account Eqs. (|25|l . (|26|) . we get standard equalities for the main correlators 

Gi"' = [t{v) ± inM"' , = gL°^ = [t^{v) + nlv^] . (34) 



I 

An explicit form of linear operator To proceed, let us consider the effective interaction 

term in action 

L = L+A + i_B + LT (35) 
obeying to equality L = (C^^^) ^ will be needed below. 

Using the equality [12 k = -T— ( ( '^^^^^^v'^v'^ 



K ^ -T- J J -^yivi\c{y^)\'\c{v^)Y (38) 
= G;^A + GI^B-G;^S'GI^T, (36) 

we arrive easily to the components 

taken in the frequency representation. Within the mean- 
Li = T(y) ± YTicV, L = — 1. (37) field approximation, one has 

I 



\c(y^)\My'i)f ^ (|c(^i)P) Wi)? + Wi)? = S(v^)\c{v^)\^ + |c(z.i)P^K) (39) 

I 

and the fluctuational component of the inhomogeneity action H38|) takes the form 

^{0} = -rA /^i/2|0(j/,,9)|2^dt? (40) 
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where parameter A given by Eq. reduces to averaged 
magnitude 



zTT / zvr 



du 



iy^C{i^,d)'0dd. (41) 



As a result, the bare mass m in the action 5*0 given by 
Eqs. 1)28(1 . H29|l is replaced by the effective quantity 



m — A 



(42) 



being averaged value of the fluctuating mass ifT^ . 

To finish supersymmetry representation of the action 
((17|l defined by the Lagrangian l|18|l . one should add to 
Eqs.|2HIl, 601) the self-action term 

V{^{z)} - / v{^{z))dz, z = {n,^} (43) 



with the expansion |(SJ). Then, the standard perturbation 
theory gives the symbolic expression p3| 

E(z?i,z?2,n) 

iP 9 d (44) 

= f|-(C^(^i,'?2;n)) +^(C(z?i,z?2;n))' 



3! 



for the self-energy function I](i?i, -(92, ") defined by the 
following equation for the n-point dressed supercorrelator 



C(") (i9, z9i)E(") (z9i , 'd2)C^°^ (t?2, t?')d^9idi?2 



(45) 



However, detailed analysis 0| shows that the multiplica- 
tion rules given by Table 1 has to be replaced by the rules 
of Table 2. Then, the components of the pseudovector 



TABLE II: 





r(i?,i?') 


A(i9,t?') 


B(i?,i9') 


r(j9,i9') 


r(i9,i?') 


A(i9,t9') 


B(i9,i9') 












B(i9,j9') 


B(i9,i9') 









S = S+A + S-B + ST 
take the following forms: 



(46) 



J 



^±(.>^) = t^'J ^5(i^i)G±(^-^i) + y 



— !— ^S'(z/i)S'(i/2)S'(j^ - i^i - 1^2), 



(2^) 
di/idi/2 



S'(i/i)5(j^2)G±(i^ - i^i - 1^2)- 



(47) 
(48) 



Making use of the theory of residues (see Appendix B) with the correlators ((34|l . where the frequency dependent 
parameter t(z^) = t — mi^^ is replaced by the bare one r, arrives at the equalities IBSp . (|B6ll which take the form 



3^ \-i 



(4r2n,)- 



A2 



A2 



A^ 



22 ^ 3^™, 

2 ;^2 



22 2 • 33rnc 



(49) 
(50) 



r 



within the hydrodynamic limit ^ = v/cos 1, lOs = t/tIc- 
Self-consistent behavior of the system under consider- 
ation is described by the generalized Dyson equation 



Combination of Eqs. jSUJ) arrives at the fi- 

nal equations for main correlators within hydrodynami- 
cal limit ^ ^ 1: 



= L S. 



(51) 



In the component representation this equality arrives at 
the equations 



5*= (E-L)G+G_, 

-1 
± 



(52) 
(53) 



T - (4T"n,)-i fi^ + 



A2 



± i 



6™, 
A2 



^-(4rV)-f^ + 



3^™ 
A2 



e (54) 



22 2-337 
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IV. DETERMINATION OF THE PERIOD OF 
MICROPHASE STRUCTURE 



that is a solution of this equation at a = /? 
result, we obtain the simple dependence 



0. As a 



Our consideration is based on the obvious equality for 
the long space period L = 21 + D where / is the oligomer 
chain length, D is the thickness of the homopolymer layer 
being fixed by the inverse share of average number of 
the hydrogen bonds (see Figure^. Physically, this value 
is reduced to the magnitude 2TT/U determined by the cir- 
cular frequency lo in the alternation of the oligomer heads 
along the homopolymer chain. Then, the long space pe- 
riod is expressed by the following equality Q (see Ap- 
pendix A) 



L = 2/ + Dow-\ Do=(2ttx 



.i/6„-i/3 



b>b (55) 



where x ^ 10^^ is the Flory parameter, n ~ 10 is the 
number of segments in oligomer chain, b is the segment 
length. 

To obtain the frequency w, one has to determine firstly 
the effective mass nief given by Eqs. H42|l . (|41|l . Usage of 
the theory of residues (see Appendix B) with the struc- 
ture factor H52|) and Green function H54|l arrives at the 
renormalization mass parameter 



A = 



1 Tlc^i^ 



1 



23 • 34 T^TUef 



(56) 



where only the terms of the second order of smallness 
over the parameters /z, A of the self-action ||SJ) are kept. 
Inserting here Eqs. ifT^ . we obtain the equations 

for determination of the effective mass as a function of 
the temperature: 



where 



T, 



cO 



wJ "2 



23-34 



r ^ ~ 25 T^n' 



(57) 
(58) 



(59) 



Numerical solution of Eq. (|58|l for different values of 
a and (3 allows us to estimate an influence of the self- 
action on the effective mass. It turned out that even 
small variation of the parameter a substantially changes 
the shape of the dependence fJ-{T), whereas the parameter 
f3 almost does not affect it, and we can put /3 = for the 
sake of simplicity. This means physically that the cubic 
anharmonicity in the self-action potential energy jSJ is 
irrelevant to the microphase separation phenomenon. 

The smallness of the self-action parameters a, (3 allows 
us to solve Eq. (|58|l analytically. In so doing, one has to 
replace the required dependence ^{T) in the right hand 
part of Eq. H58|) by the bare dependence 



A'o(T) 




(60) 



M(T) 



1 



(61) 



(62) 



with a characteristic temperature 

= Teo (1 - 2a) 



where the scale Tco is given by the first of Eqs. l|5^ 
(the multiplier should be put /io(7co) ~ 1/2 due to the 
smallness of the parameter a ^ 1). According to Eq. 
(|61|l . with the temperature increasing the effective mass 
(|57|l decreases monotonously from the bare magnitude m 
at T = to m/2 at r = Tc (see the main panel in Figure 
121. The critical temperature Tc determines the point of 
the order-disorder transition according to the condition 

d/i 

-TT^ = —CO- 
AT 

Resulting dependence Tc on the self-action parameter a 
is shown in the inset of Figure |21 It is principally impor- 
tant that the bigger value of the self-action parameter a, 
the more narrow is the temperature domain Tg-^Tc where 
the microphase separated structure is possible {Tg being 
a glassing temperature). In other words, the self-action 
effect leads to the shrinking the region of the ordered 
structure because the critical temperature Tc reaches the 
boundary magnitude Tg with increasing of a before the 
magnitude a ~ 0.08. 




FIG. 2: Temperature dependence of the inhomogeneity pa- 
rameter /i for different values of the self-action parameter 
a: thick curve relates to a = 0, thin one — to a = 0.08 
(Tco = 337K). Inset: the temperature of order-disorder tran- 
sition Tc versus the parameter a. 

The divergency condition GZ^ = of the Green func- 
tion (|54|) gives the proper frequency 



1^0 = — ic!7, uj — y uJq — vo'^ (63) 
of the oligomer alternation along the homopolymer chain. 
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Real and imaginary parts are determined by the expres- 
sions 



Wo 



Wo 



<^o To 
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KT) T 







1 + 6a (^1 + ^ 






n). 



(64) 



where the dependence ^(T) is defined by Eqs. l|HT|) . 

the effective mass in parenthesis after the factor a ^ 1 



is put to be equal to the value m/2 related to the critical 
temperature T^; characteristic scales of both frequency 
and temperature are introduced as follows: 




(65) 



As a result, combination of Eas. (|55|l . (|63|l and H64|) leads 
to the final result for the long space period 



3 1 + 8(To/Tf 
2" 1 - 2(To/T)2 



(To/T)- 



(66) 



r 



where the characteristic length Lq — Dq/luq ^ 
y^fh/r x^^^b oc x^^^N^^'^ is the function of both param- 
eters X Etnd A'' being thermodynamically independent. 
Thus, the first of the exponents in the scaling relation 
Lq oc x"7V'' takes the magnitude a = 1/6 inherent to 
the strong segregation regime, whereas the second one 
{b ~ 1/2) is the same for the weak one 19]. Note that 
the obtained x-dependence is caused by the multiplier 



X 



1/6 



in the generic relation (|55|l that is relevant to the 
former of above regimes, while the method developed ad- 
dresses to the latter. 



V. DISCUSSIONS 

The behavior of the system under consideration is con- 
trolled by the parameters m, r and a which determine the 
temperature Tc of the order-disorder transition and the 
long space period L given by Eqs. (|^ . respectively. 
Moreover, there is the self-action parameter < a <C 1 
whose value is limited by the magnitude amax — 0.08 
(see insert in Figure |2J). To guarantee positive values 
of the radicand in Eq. H66I) at the critical temperature 
Tc, the above parameters have to be constrained by the 
condition 

K>V2 

limiting magnitudes of the principal parameter 

Tc 

K — — - 

To 

The minimal magnitude of k fixes the choice of the theory 
parameters according to the condition 

(1 



2a). 



(67) 



(68) 



a < 2-i/4m5/V/^ 



(69) 



It would seem from Eqs. (|67|l . (|68|l that the decrease 
of the critical temperature Tc with passing from the ioni- 
cally bonded system (such as P4VP-(DBSA)j.) to the hy- 
drogen bonded one (e.g., P4VP-(PDP)2.) is caused only 



by the growth of the fluctuation parameter a with re- 
spect to the mean magnitude of the inhomogeneity pa- 
rameter ffi. It appears, however, that the main reason 
for such behavior is given by the decrease of the mean- 
geometrical magnitude mr of the principal coefficients 
in the generic Lagrangian © (see below). 




FIG. 3: Long space period in the strongly bonded systems as a 
function of the oligomeric fraction x. Solid lines represents the 
results of fitting in accordance with Eq. Experimental 
data for P4VP-(DBSA)^ (•) and P4VP-(Zn(DBS)2)^ (o) at 
room temperature are taken from Ref. 0- 

To clarify this problem and find explicit form of the 
dependencies of the temperature of order-disorder tran- 
sition Tc and the period L on the oligomeric fraction 
X, we assume for main theory parameters m and r the 
three-parametric relations: 

m = mo + Ax{xm — x), t = tq + Bx{xr ~ x) (70) 

with positive constants too, 'Ha, ^, 5, Xr to be de- 
termined. Then, the fitting of the experimental results 
shown in the Figure |3| in accordance with Eq. H66|) where 
T, TO are given by Eq. (|70|l leads to the following results 
for the ionically bonded systems: 
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FIG. 4: Order-disorder transition temperature Tc for the 
weakly bonded system as a function of the oligomeric frac- 
tion X. The sohd hne represents the dependence obtained by 
fitting according to Eq. 1621 1. Experimental data for P4VP- 
(PDP)^ (■) are taken from Ref. Q. 



FIG. 6: Temperature dependence of the long space period 
in the weakly bonded system. The solid line represents the 
dependence obtained by fitting according to Eq. Ifcibl . Ex- 
perimental data for P4VP-(PDP)a; at x — 0.85 (o) are taken 
from Ref. |||. 




0,0 0,3 0,6 0,9 1,2 1,5 1,8 2,1 

X 



FIG. 5: Long space period in the weakly bonded systems as 
a function of the oligomeric fraction x. Solid lines represents 
the results of fitting in accordance with Eq. I|66^ . Experimen- 
tal data for P4VP-(PDP)a; at temperature of order-disorder 
transition Tc (■) and at temperature T — 80°C (□) are taken 
from Ref. ||]. 

• the mixtures P4VP-(DBSA)2; 

mo = 18, A — 8, Xm — 1.5; 

To = 0.6, B = 1.5, Xr^l-O; (71) 

a — 0.01; b — Inm; I — lOnm; 

• the mixtures P4VP-(Zn(DBS)2)a; 

mo = 5.3, A = 26, x„i = 1.6; 

To = 0.8, B = 0.1, Xr = 1.0; (72) 

a = 0.01; b = Inm; / = lOnm. 

At X = 1 one obtains m = 22, r = 0.6 for P4VP- 
(DBSA)^; and m = 20.9, r = 0.8 for P4VP-(Zn(DBS)2):r. 



Then, Eq. ^ gives values k = 10^, 10^ at cr = 1.31, 
4.18, respectively. 

Much more complicated situation occurs in the weakly 
bonded system P4VP-(PDP)a;. Here, decrease of the pa- 
rameter (|68f) results in the narrowing of the temperature 
domain Tq Tc of the phase separation. All parame- 
ters for this class of systems can be determined by the 
combined fitting of a series of experimental data for the 
critical temperature Tc and the long space period L (see 
Figures 0] — First constraints follow from the com- 
parison of experimental points for the temperature Tc of 
order-disorder transition (see Figure ^ with fitting re- 
sults based on Eq. at a = 0.01, I = 10 nm: 

^r^ = 562, — = 0.155, x„ = 1.615. 73) 

cr^ Too 

The following of parameters gives application of Eq. H66(l 
for the long space period at the temperature T — Tc to 
the data shown in Figure|31as the non-monotonous curve: 



— = 1499, — = 7.968, Xr = 1.926. (74) 

To To 

Finally, making use of the expression (|66|l and experi- 
mental data for the temperature dependence of the long 
space period given in Figure yields the last constraint 

, ^ = 924. (75) 

As a result, taking too = 1 at x = 1 the magnitudes 
A = 0.155, e = 23.87, a = 0.206 are obtained to provide 
extremely small value r = 5.6-10"'^ of the hydrogen 
bonding strength and the temperature scale To = 160. At 
a = 0.01 this arrives to the rest of parameters k = 2.07, 
A = 6.22- 10-^. 



10 



TABLE III: 





ST '-t \ IT I J_^J_> O-TA- J X 


P4VP— ("/n /'DRS^-.'! 


P4VP— fPDP^ 

IT ^ \ ST \ ST i—f iT J 


TTlQ 


18 


5.3 


1.000 


A 


8 


26 


0.155 


■^m 


1.5 


1.6 


1.615 


TTi 


22 


20.9 


1.095 




0.6 


0.8 


6.67 ■ 10~" 


B 


1.5 


0.1 


5.3 ■ 10~^ 


Xt 


1.0 


1.0 


1.926 


f 


0.6 


0.8 


5.6 ■ 10"^ 


(J 


1.31 


4.18 


0.206 


OL 


0.01 


0.01 


0.01 


X 


3.33 


6.66 


6.22 ■ 10"" 


e,K 


10'' 


lO'' 


23.87 




1.41 • lO'' 


1.25 ■ 10" 


337 




1.38 • lO'' 


1.22 ■ 10" 


331 


To,K 


138 


122 


160 


K 


10^ 


10^ 


2.07 


I, nm 


10 


10 


10 


b, nm 


1 


1 


1 



It is worthwhile to discuss separately the dependence 
of the long space period on the oligomer/monomer ra- 
tio at the temperature T = 80" C that relates to the 
monotonous decaying curve shown in Figure IS] Because 
the maximal temperature of the order-disorder transi- 
tion is Tc ~ 65°C to be corresponded to a; = 0.85 (see 
Figure 0J), experimental data related to T = 80°C are 
obtained for the temperature being beyond of the region 
of the ordered structure (T > Tc). From the physical 
point of view, at the critical temperature T = Tc the pe- 
riodicity of the microphase separated structure formed is 
caused by long-range correlations, whereas at T = 80"C 
only short-range correlations hold to be determined by 
the homopolymer backbone together with the hydrogen 
bound surfactant molecules P| . Fitting of the experimen- 
tal points for the dependence L{x) at the temperature 
T = 80^0 can be done on the base of Eq. (|66|l where one 
puts /i(T) — fJ.{Tc) = 1/2. Then, the values of the pa- 
rameters obtained differ from those obtained for T = Tc 
by the following constraints: 

B 

_ = 4.2, xr = 3.27, , = 543. (76) 

To ^JmoTo 

Obviously, this difference is due to the temperature de- 
pendence of the hydrogen bonding parameter r in the 
potential energy 

To conclude our estimations, we notice the model 
developed explains successfully a vast variety of pecu- 
liarities obtained experimentally for various classes of 
homopolymer-oligomer mixtures with the interactions of 
different strength. The resulted data for strong, interme- 
diate and weak coupled systems P4VP-(DBSA)a;, P4VP- 
(Zn(DBS)2)^ and P4VP-(PDP)^, respectively, are given 
in Table III. It is seen that the coupling weakening ar- 
rives to a decrease of both inhomogeneity parameters m 



and cr, as well as to the crucial decrease of the hydro- 
gen bonding parameter r and the self-action parameter 
A, on the one hand, and the characteristic temperatures 
Tc and 8, on the other hand. According to the relations 
(|68|1 this leads to extremely large suppression of the value 
of the parameter k, that causes the crucial shrinking the 
temperature interval of the microphase separation. An 
analogous effect is caused by the self-action increase. 

To get rid of a misunderstanding, we would like to 
stress a composite character of the approach used. As 
it is mentioned in Introduction, this circumstance is ex- 
pressed by dividing the total free energy ^ into two 
terms, the first one Fho is relevant to the non-associated 
homopolymer-oligomer mixture, the second one Fhb is 
caused by the hydrogen bonding. These terms are caused 
by the interactions of principally different physical na- 
ture: the behavior of the mixture of non-associated ho- 
mopolymers and oligomers is determined by the Flory 
parameter characterizing unfavorable interactions be- 
tween the the oligomers and the rest of the system; the 
temperature induced distribution of hydrogen bonds is 
determined by the parameter r, giving the strength of 
this bonding. From the formal point of view, both of 
the above contributions Fho(Xi4') and Fhb(j,x) should 
have a similar dependencies on the state parameters be- 
ing (apart from the temperature) the volume fraction of 
the homopolymer (/> for the first contribution, and the 
oligomer/monomer ratio x for the second one. Because 
the term Fho ~ x4'{'^ ~ 4') involves the parabolic depen- 
dence on the parameter bounded by maximal value 
(j) — 1, we took generalized parabolic approximation l|70|) 
for the dependence of the hydrogen bonding strength r 
on the oligomer/monomer ratio x which may take values 
X > 1. 

Apart from the above difference in nature of the inter- 
actions, one needs to emphasize at once the difference in 
the approaches used: the mixture of non-associated ho- 
mopolymers and oligomers had been studied within the 
strong segregation limit 5j j whereas for the consideration 
of the hydrogen bonding we use opposite approach. This 
difference is kept if the Flory parameter takes large val- 
ues X 5: 10"^, whereas the hydrogen bonding strength is 
relatively small (r <C 1). Indeed, the formula l|55() for the 
long space period was obtained within approximation of 
the sharp interface, which thickness is A ~ x~^^'^b > 3b 
to be relevant to the strong segregation regime (5|. In 
the consideration presented, we have focused mainly on 
the study of the hydrogen bonding on the base of the 
action ^ that has the form of series in powers of the 
order parameter c and its derivatives c. Such an expan- 
sion supposes making use of the weak segregation limit 
corresponding to the small values of the parameters m 
and T. 

Finally, it is worthwhile to discuss a difference with an 
usual picture of the phase transitions that is caused by 
the self-consistency condition H41|) . A critical value of 
the Flory parameter Xc in usual copolymers is known to 
be caused by the self-action effects. The accounting of 
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these effects arrives to replacement the bare parameter x 
by the renormalized value X~ Xc However, in our 

case the value of Flory parameter is so large that the tem- 
perature of the separation of non-associated polymer- 
oligomer mixture is negligible small. As a result, the 
role of X passes to the hydrogen bonding parameter r 
which does not relate to the tendency of monomers of 
the different kinds to avoid each other. However, as it is 
shown by the considerations given in , , understand- 
ing of the whole picture of microphase separation, includ- 
ing the temperature dependence of the structure period, 
demands accounting the inhomogeneity in the distribu- 
tion of oligomers along homopolymer chains. Within the 
approach developed, this is reached by means of the effec- 
tive kinetic energy with the mass fluctuating due to 
the temperature dependence of hydrogen bonding. This 
dependence leads to the reduction (|42() of the effective 
mass mrief that causes a phase transition from stochastic 
to periodic distribution of the oligomers along the ho- 
mopolymer chain. However, if the critical point is fixed 
usually by the condition nief = [T]| . in our case the 
critical temperature Tc relates to the finite magnitude 
nief = m/2 of the effective mass which has a singular- 
ity dme//dT = — oo in the temperature derivative (see 
Figure 2). 
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APPENDIX A: DERIVATION OF A GENERIC 
RELATION FOR MICROPHASE STRUCTURE 
PERIOD 

Following we suppose the period of the microphase 
structure to be determined by the minimum of the spe- 
cific free energy 



1 Pint + Fs; 



V 



T 



V = LS, L=2l + D (Al) 



related to the first term in Eq. . Being the free energy 
of the homopolymer-oligomer mixture, this term consists 
of the interfacial and stretching components i^mt, -Fstr 
measured in the temperature units T per the domain 
volume V (according to Figure 1 L, I and D are the 
long space period, the length of the oligomer tail and the 
thickness of the homopolymer layer, respectively; S is the 
domain surface area). 

The interfacial free energy is stipulated by the loss of 
conformational entropy caused by the localization of the 
homopolymer chains within the interface of thickness A. 
Due to unfavorable interaction x between the oligomer 
tails and the polymer layer the chains form up loops con- 
taining segments of number A/'s ~ x~^ E3- Then, within 



the model of the random walk, the interface thickness is 
estimated by the relations ~ A/'s 6^ ~ b'^/x where 
b is the segment length. Respectively, the interfacial 
free energy Fint — Mi T is determined by the number 
Afi ~ SA/Afgb^ of the loops within the interface. As a 
result, we obtain the estimation 



X 



1/2 



52 



-T S. 



(A2) 



Another addition F^tr is caused by the stretching of 
the surfactant side chains, whereas the stretching of 
the homopolymer chains enlarges only the volume part 
of the free energy. This addition is expressed by the 
simple equality F^tr — J^c^sFi where the first factor 
Afc ~ DS/Nb^ gives the number of chains per layer, 
the second multiplier Hs ~ (b/X)N is the number of the 
oligomer molecules per chain of N segments (A is period 
of the oligomers alternating along the chain) and the last 
factor Fi ~ {l'^/nb'^)T presents the free energy of stretch- 
ing a side chain of n segments to a length I. Combining 
the above multipliers, we find the estimation for the total 
free energy of stretching 

F„~^T^ (A3) 

To derive the explicit expression for the dependence of 
the free energy l|Aip on the layer thickness D we need to 
use an obvious condition 21S = JKfcUsVs where Vs — nb^ is 
the volume of surfactant molecule. As a result, we obtain 
the relation 



21 b 



(A4) 



according to which the period A defines the rest of 
geometrical characteristics of the microphase separated 
structure. Inserting Eqs. ljA2|l — (|A4|I into Eq. ljAl|l . 
we arrive to the final expression for the interfacial free 
energy: 



/ 



(A5) 



[1 + {b/X)n] 62 

where numerical coefficients are dropped. The minimiza- 
tion condition df /dD = arrives at the steady-state val- 
ues of the homopolymer layer thickness and the oligomer 
length 

.1/6 



D 



X 



,1/3 



A, 21 - b. (A6) 

It is convenient to express above results by means of 
the dimensionless frequency of the oligomer alternating 
along the homopolymer chain: 



2tt 
Xjb 



r.X- 



27r, A = - (A7) 
A 



where X is the averaged oligomeric fraction per ho- 
mopolymer. So, the long space period of the microphase 
separated structure takes the form 



L = 2l + D'. 



(A8) 
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APPENDIX B: CALCULATIONS OF 
CONVOLUTION INTEGRALS 

a. Self-energy functions Calculations of the self- 
energy functions l|T7|) . lead to a rather tedious pro- 
cedure due to the convolution integrals. To demonstrate 
the line of these calculations we consider in details the 
simplest integral related to the first term in Eq. H48|) : 



(Bl) 



Making use of the expressions H34|l for the structure fac- 
tor S and Green function G-, where the frequency dis- 
persed parameter t{i/) is replaced by its bare magnitude 
r, arrives at the convolution integral 



(B2) 



where a characteristic frequency Ug = r/ric is introduced. 
This integral has the poles iiw^ and v 4- icj^ (see Figure 
7a). In accordance with the theory of residues the 











ic, 






a) -''^s 


V -6 
b) 
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expression yields 



V 



(B4) 



Analogously, the rest of convolution integrals is calcu- 
lated giving final form of Eqs.(|37I), (|^ : 



E = 



1 



A2 



1 



E+ - 



2t3; 



4 + ^2 

4 + ^2 



8r4n2 9-1-^2' 

3 Tig 
8t3?i2 9 + ^2 ■ 



(B5) 
(B6) 



&. Renormalization mass parameter Explicit form 
of the renormalization mass parameter H41|l is determined 
by the structure factor (|52|l and Green function H54|l with 
the effective mass (I42II and parameter t{i/) being replaced 
by bare r: 



A = 



■KTlcmef J (^2 

where one denotes 



2ieie-el)(e2 + 2iae-el 



(B7) 



^0 

6 



IGrriefT 

He 
2mefT 

n2 



1 



1 



A2 



4r3r 



9tti, 
A2 



(B8) 



FIG. 7: Poles of the convolution integral in Eqs. IIB2II (a) and 



integral in Eq. (|B2|) is reduced to sum over two of these 
residues that locate in upper half-plane of the complex 
frequency vi. 



27ri 



2iws • {—v) V ■ (y + 2\bJs) 



(B3) 



where terms in the square brackets relate to the poles 
iujg and v -f- iLOg, respectively. After a simple algebra this 



The integral in Eq. (|B7() has the pole structure that is 
shown in Figure 7b. As above, the sum over residues 
located in the upper half-plane of the complex frequency 
^ yields the integral value 



2 1 Coa 



(B9) 



With accounting the notices l|B8l) and keeping only the 
terms of the second order of smallncss over parameters fi 
and A, one obtains the final expression 15t)|l . 
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